A Parameterized Post-Priedmann Framework for Modified Gravity 
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We develop a parameterized post-Friedmann (PPF) framework which describes three regimes of 
modified gravity models that accelerate the expansion without dark energy. On large scales, the 
evolution of scalar metric and density perturbations must be compatible with the expansion history 
defined by distance measures. On intermediate scales in the linear regime, they form a scalar-tensor 
theory with a modified Poisson equation. On small scales in dark matter halos such as our own 
galaxy, modifications must be suppressed in order to satisfy stringent local tests of general relativity. 
We describe these regimes with three free functions and two parameters: the relationship between 
the two metric fluctuations, the large and intermediate scale relationships to density fluctuations 
and the two scales of the transitions between the regimes. We also clarify the formal equivalence 
of modified gravity and generalized dark energy. The PPF description of linear fluctuation in f(R) 
modified action and the Dvali-Gabadadze-Porrati braneworld models show excellent agreement with 
explicit calculations. Lacking cosmological simulations of these models, our non-linear halo-model 
description remains an ansatz but one that enables well-motivated consistency tests of general 
relativity. The required suppression of modifications within dark matter halos suggests that the 
linear and weakly non-linear regimes are better suited for making complementary test of general 
relativity than the deeply non-linear regime. 

PACS numbers: 



I. INTRODUCTION 

Theoretically compelling alternatives to a cosmological 
constant as the source of the observed cosmic acceleration 
are currently lacking. In the absence of such alternatives, 
it is useful to have a phenomcnological parameterized ap- 
proach for testing the predictions of a cosmological con- 
stant and phrasing constraints in a model-independent 
language. This approach parallels that of local tests of 
general relativity. The parameterized post-Newtonian 
description of gravity forms a complete description of 
leading order deviations from general relativity locally 
under a well-defined set of assumptions 

A parameterization of cosmic acceleration from the 
standpoint of dark energy is now well-established. The 
expansion history that controls distance observables is 
completely determined by the current dark energy den- 
sity and its equation of state as a function of redshift. 
Structure formation tests involve additional parameters 
that control inhomogencities in the dark energy. Covari- 
ant conservation of energy-momentum requires that the 
dark energy respond to metric or gravitational poten- 
tial fluctuations at least on scales above the horizon. In 
a wide class of models where the dark energy remains 
smooth relative to the matter on small scales, the phe- 
nomenological parameter of interest is where this transi- 
tion occurs @, 13, 0| • 

A similar structure is imposed on modified gravity 
models that accelerate the expansion without dark en- 
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ergy. Requirements that gravity remain a metric theory 
where energy-momentum is covariantly conserved also 
place strong constraints their scalar degrees of freedom. 
On scales above the horizon, structure evolution must be 
compatible with the background expansion [5j. Interme- 
diate scales are characterized by a scalar-tensor theory 
with a modified Poisson equation ((J. If these modifi- 
cations are to pass stringent local tests of gravity then 
additional scalar degrees of freedom must be suppressed 
locally 7]. Two explicit models that exhibit all three 
regimes of modified gravity are the so-calledf (R) mod- 
ified Einstein-Hilbert action models [g, 0, Qjj] and the 
Dvali-Gabadadze-Porrati (DGP) braneworld model [TTI | . 

Although several parameterized gravity approaches ex- 
ist in the literature, none describe all three regimes of 
modified gravity (cf. [2 QS 03) 

and most do not ex- 
plicitly enforce a metric structure to gravity or energy 
momentum conservation (e.g. [IM [Trj. \vj\, [Lsl. Il9j ) . 

In this paper, we develop a parameterized post- 
Friedmann (PPF) framework that describes all three 
regimes of modified gravity models that accelerate the 
expansion without dark energy. We begin in [jlT] by de- 
scribing the three regimes individually and the require- 
ments they impose on the structure of such modifications. 
In ^IIIl we describe a linear theory parameterization of 
the first two regimes and test it against explicit calcula- 
tions of the f(R) and DGP models. In mVl we develop a 
non-linear ansatz for the third regime based on the halo 
model of non-linear clustering. In the Appendix, we clar- 
ify the formal relationship between modified gravity and 
dark energy beyond the smooth class of models. 
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II. THREE REGIMES OF MODIFIED GRAVITY 

In this section, we discuss the three regimes of modified 
gravity theories that accelerate the expansion without 
dark energy. We begin by reviewing the requirements on 
super-horizon metric perturbations imposed by compat- 
ibility with a given expansion history ( W Ap . Modifica- 
tions that introduce extra scalar degrees of freedom then 
enter a quasi-static regime characterized by a modified 
Poisson equation ( 311 Bp . Finally stringent local tests of 
gravity require that modifications are suppressed in the 
non-linear regime within collapsed dark matter haloes 

(3E3- 



A. Post-Friedmann Super-horizon Regime 

Under the assumption that modified gravity remains 
a metric theory in a statistically homogeneous and 
isotropic cosmology where energy-momentum is covari- 
antly conserved, a parameterization of the expansion his- 
tory that is complete under general relativity is complete 
under modified gravity as well. A modified gravity model 
and a dark energy model with the same expansion rate 
H = a~ 1 da/dt and spatial curvature predicts the same 
observables for any measure that is based on the distance- 
redshift relation (see e.g. Hence in terms of the 

background, modified gravity models can be parameter- 
ized in the same way as dark energy without loss of gen- 
erality. Neglecting spatial curvature and radiation for 
simplicity here and throughout, we can assign an effec- 
tive energy density 
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would be the contribution of matter to the expansion 
under the normal Friedmann equation. Alternately we 
can assign a current effective energy density in units of 
the critical density f2 ff = 8ttGp c s/3Hq and an effective 
equation of state 
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Compatibility with this expansion combined with en- 
ergy momentum conservation highly constrains the evo- 
lution of metric fluctuations above the horizon. Super- 
horizon metric fluctuations in a perturbed universe can 
be viewed as evolving as a separate universe under the 
same modified Friedmann equation but with different pa- 
rameters. 

Bertschinger [f| showed that consequently metric fluc- 
tuations in fact obey the same fundamental constraints 



as they do in general relativity. These constraints ap- 
pear in different ways in different gauges as detailed in 
the Appendix. Under the assumptions of a metric the- 
ory and energy-momentum conservation, all of the usual 
gauge structure including the so-called "gauge invariant" 
approach used here apply to modified gravity as well. 

In the comoving gauge of the matter, the constraint 
for adiabatic initial conditions looks particularly simple. 
The curvature or space-space piece of the metric fluctu- 
ation £ remains constant to leading order ([2lj], see also 



C = 0{k 2 H Q 



(4) 



where ' = d/dhia and ku = k/aH is the wavenumber 
in units of the Hubble parameter. In the more famil- 
iar Newtonian gauge where the curvature is denoted $ 
and the time-time piece or gravitational potential 'J, the 
gauge transformation equation (|A19|) 



C = $ - V m /k H 
and the momentum conservation equation 

+ V m = k H ^> 
along with Eqn. ([4]) imply 
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Here V m is the scalar velocity fluctuation of the matter in 
both the comoving and Newtonian gauges (see Eqn. (|A5|) 
and Eqn. (IA10|) V Eqn. ([7]) is also satisfied in general 
relativity 

These relations have been explicitly shown to hold for 
DGP braneworld gravity [23j and f(R) modified action 
gravity [24]]. What distinguishes a particular model of 
gravity or dark energy is the relationship between <E> and 
^ in the Newtonian gauge or equivalently £ and V m in 
the comoving gauge. We will parameterize this relation 
by 
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In terms of the post-Newtonian parameter 7 = — ^/vP, 
g = (7 — l)/(7 + 1). Given that gravitational rcdshift 
and lensing effects involve the metric combination 
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we will typically state metric results in terms of It 
is useful to note that $ = (g + 1)$_ and 9 = (g — 1)<3>_. 
Super-horizon scalar metric fluctuations for adiabatic 
perturbations are completely defined by the expansion 
history H and the metric ratio g. 

In Fig. [T] we show the evolution of given a metric ra- 
tio that evolves as g = g$a for a ACDM expansion history 
defined by u> G ff = —1 and 1 — f2 c ff = f2 m = 0.24. Given 
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FIG. 1: Super-horizon (SH) metric evolution of <E>_ = ($ — 
^)/2 for various choices of g = goa. For g < 0, 5>_ can 
actually grow near the onset of acceleration. The expansion 
history is fixed by to e ff = — 1 and Q m = 0.24 in all cases. 



that g — > as a — *■ 0, wc take matter dominated initial 
conditions of <&_ = — 3£i/5 and = as in gen- 
eral relativity. Note that if g < 0, <E>_ can actually grow 
during the acceleration epoch. Changing the evolution of 
gravitational potentials alters the low order multipoles of 
the CMB through the integrated Sachs Wolfe effect. For 
example it can be used to suppress the CMB quadrupole 
(see e.g. (H). 

Under general relativity, the metric ratio is determined 
by the ratio of anisotropic stress to energy density (see 
Eqn. (|A2ip ). It supplements the information in the ex- 
pansion history since in the background the anisotropic 
stress vanishes by assumption. In particular, with the 
matter carrying negligible anisotropic stress, this rela- 
tionship is determined by the anisotropic stress of the 
dark energy. For models of dark energy based on scalar 
fields, this also vanishes in linear theory and hence g = 0. 
Nevertheless, a hypothetical dark energy component that 
produces the same g as a modified gravity model cannot 
be distinguished from super-horizon metric fluctuations. 
We further examine the relationship between modified 
gravity and dark energy in the Appendix. 



B. Post-Newtonian Quasi-Static Regime 

Well inside the horizon but still in the regime of linear 
fluctuations, theories that exhibit an extra scalar degree 
of freedom tend toward a post-Newtonian scalar-tensor 
description. As is the case in general relativity when time 
derivatives of the metric fluctuations can be ignored com- 
pared with spatial gradients, the modified field equations 
reduce to a modified Poisson equation 
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where A m is the fractional density perturbation (in the 
matter comoving gauge), fa parameterizes a possibly 
time-dependent modification of the Newton constant, 
and the relationship between the two metric fluctua- 
tions is again parameterized by g as in Eqn. ([8]). We 
assume that spatial fluctuations in Jq lead to a second 
order correction in the linear regime. We call this the 
quasi-static approximation. This quasi-static approxi- 
mation has been explicitly shown to hold for both DGP 
braneworld gravity [23|, [25| and f(R) models [HIHH once 
k H > 1. 

Density fluctuations on the other hand are determined 
by the Newtonian limit of the conservation equations 
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- V m = k H V = {g - l)k H $- 
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Combining the modified Poisson equation (|T0|) and the 
conservation equations and taking Jq =const. for 
simplicity yields 



~W 

H' 3 
" ~H + 2 (T 



f G )H*d> 



-(.9-1) 



(12) 
= 0. 



Note that this quasi-static equation is inequivalent to the 
super-horizon evolution ([7]) whenever w c s ^ — 1, g ^ 0, 
or /g 7^ 0. The differences for w c g ^ — 1 applies in gen- 
eral relativity as well and corresponds to the well-known 
fact that growth defined by a smooth dark energy com- 
ponent is inconsistent with conservation of energy mo- 
mentum on super-horizon scales (e.g. @, Q). In Fig. [2 
we compare the quasi-static and super-horizon evolution 
to the present. Note that in terms of the change in the 
gravitational potential from its initial value, important 
for gravitational redshift effects in the CMB, the dif- 
ferences are amplified. Compared with ACDM where 
A$_ i=s — $i/4 (see Fig. [TJ, these changes are enhanced 
by a factor of ~ 4 and cannot be neglected for large values 
of| 5 |. 



C. General Relativistic Non-linear Regime 

A successful modification of gravity must have a 
third regime where non-linearities in the modified field 
equations bring the dynamic back to general relativ- 
ity. Scalar-tensor modifications of gravity that persist 
to small scales can be ruled out by stringent local tests 
of gravity. For example, the Cassini mission imposes the 
limit Q 



\g\ < 1.2 x 10" 



(13) 



for the metric ratio in the solar system. 

In both the DGP braneworld and f(R) modifications, 
the extra scalar degree of freedom obeys a non-linear 
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FIG. 2: Fractional difference between the quasi-static (QS) 
vs super-horizon (SH) metric at a = 1 for a metric ratio of 
g = goa. The two evolutions differ unless w c ff = — 1 and 
go — 0. Q m — 0.24 for all cases. 



equation that suppresses its effects within collapsed ob- 
jects such as dark matter halos. In the DGP braneworld 
model, the scalar degree of freedom represents a brane 
bending mode. Non-linear interactions of this mode be- 
come important at the so-called Vainshtein radius r* ~ 
(2GMr^) 1 / 3 Q- Here r c ~ H^ 1 is the crossover scale 
where gravity becomes fully 5D. This radius is compara- 
ble to the virial radius of dark matter halos and suggests 
that inside a halo gravitational interactions behave as in 
general relativity. 

Similarly, for modified action f(R) models the extra 
scalar degree of freedom corresponds to df /dR and has a 
mass that depends on the local curvature R. Deviations 
from general relativity can then be sup pressed by the 
so-called chameleon mechanism [2?], HI] so long as the 
gravitational potential is sufficiently deep [1^, H3| . 

In the empirical PPF description that follows in ^IIII 
and ^IVI we seek a description that joins these three 
regimes. 



III. LINEAR THEORY PARAMETERIZATION 

In this section, we construct a parameterized frame- 
work for linear perturbations in modified gravity models 
that joins the super-horizon and quasi-static regimes de- 
scribed in the previous section. We first describe the 
construction ( ^III Al) and then test the parameterization 
against modified action f(R) models fi jlll B[) and the 
DGP braneworld model ( fllTTC) . 



A. PPF Parameters 

We have seen in fJTT] that in both the super-horizon and 
quasi-static regimes, the evolution of metric fluctuations 
are primarily determined by the metric ratio g. How- 
ever, the manner in which the metric ratio g determines 
metric evolution differs between the two regimes. The 
quasi-static regime also allows the freedom to change the 
effective Newton constant. 

Let us introduce a parameterization that bridges the 
dynamics of the two regimes at a scale that is parameter- 
ized in units of the Hubble scale. The Newtonian-limit 
conservation equations (jlip must first be corrected for 
metric evolution. The exact conservation equations im- 
posed by V^Tf+u = and the metric is given by (see 
Eqn. (jA15|) and (|A20l) ) 



A' m = -k H V m - 3C , 
VL + V m = (g- l)fc H $_ , 



(14) 



where the additional term involving evolution of the met- 
ric is given in the matter comoving gauge £' to match the 
definition of density perturbations A m in this gauge. It 
is related to the evolution of the Newtonian metric by 
Eqn. © 

C = (g + l)&_ + (l-g + g')$--^^. (15) 

H k H 

In order to match the super-horizon scale behavior we 
introduce an additional term T to the modified Poisson 
equation dTUJ) 



/c 2 [$_ +T] = AnGa 2 p m A r 



(16) 



We now demand that as kn — ► T enforces the met- 
ric evolution of Eqn. ([?])• In this limit, the derivative 
of Eqn. (|16p gives an evolution equation for T given the 
conservation equations f| L4|) and the required metric evo- 
lution. Aside from g, the only remaining freedom is de- 
termining the leading order behavior of £' . Without loss 
of generality, we can parameterize Eqn. (U) with a possi- 
bly time-dependent function /<; 
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Although the super-horizon metric is determined by H 
and g alone, its relationship to the comoving density 
perturbation is not. Since k H V m = ©(fc^C) an d C = 
0(^0, this degree of freedom enters into the conser- 
vation equation (| L4[) at leading order. Combining these 
relations, we obtain the equation of motion for T 



r' + r = s, (k H ^o), 

where the source is 

1 H 3 H vn /1 » \ 
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Here we have kept only the leading order term in Uh- 

Note that the exact choice of /<; is rarely important 
for observable quantities. Any choice will produce the 
correct behavior of the metric evolution since that de- 
pends only on enforcing £' = 0(k 2 1 C,). Hence observ- 
ables associated with gravitational redshifts and lensing 
are not sensitive to this choice. Only observables that 
depend on the comoving density on large scales beyond 
the quasi-static regime are affected by this parameter. 
Furthermore the super-horizon density perturbation in 
Newtonian gauge or any gauge where the density fluctu- 
ation evolves as the metric fluctuation is also insensitive 
to f c . 

On small scales, recovery of the modified Poisson equa- 
tion (fit)]) from (TTT))) implies 

r = /c*- , (k H -» oo) . (20) 

Finally to interpolate between these two limits we take 
the full equation of motion for T to be 

(1 + 4k%) [T' + T + c\k\ (r - / G $_)] = S . (21) 

For models where S — > as a — > we take initial con- 
ditions of r = r' = when the mode was above the 
horizon. 

In summary, given an expansion history H{a) 1 our 
PPF parameterization is defined by 3 functions and 1 pa- 
rameter: the metric ratio <7(lna,fci/), the super-horizon 
relationship between the metric and density /^(lna), the 
quasi-static relationship or scaling of Newton constant 
/c(lna), and the relationship between the transition 
scale and the Hubble scale cr- For models which modify 
gravity only well after matter radiation equality, these 
relations for the metric, density and velocity evolution 
combined with the usual transfer functions completely 
specify the linear observables of the model. In specific 
models, these functions can themselves be simply param- 
eterized as we shall now show for the f(R) and DGP 
models. 



B. f(R) Models 
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FIG. 3: Evolution and scale dependence of the metric ratio g 
in f(R) models compared with the PPF fit. Here Bo = 0.4, 
w c g = — 1 and Q m = 0.24. 
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FIG. 4: Evolution and scale dependence of $_ in f(R) models 
compared with the PPF fit. Here Bo = 0.4, w e ff = — 1 and 
On = 0.24. 



In f(R) models, the Einstein- Hilbcrt action is supple- 
mented by the addition of a free function of the Ricci 
scalar R. The critical property of these models is the ex- 
istence of an extra scalar degree of freedom fn = df/ dR 
and the inverse-mass or Compton scale associated with 
it. The square of this length in units of the Hubble length 
is proportional to 



B = 



fRR p/ H 

1 + fn H' 



(22) 



where fuR = d 2 f/dR 2 . Below the Compton scale, the 
metric ratio g — > —1/3. 

The evolution of B and the expansion history come 
from solving the modified Friedmann equation obtained 



by varying the action with respect to the metric. We 
follow the parameterized approach of [24[ where a choice 
of the expansion history through w e ff and the Compton 
scale today Bq = B(\na = 0) implicitly describes the 
f(R) function and model. For illustrative purposes, we 
take f2 m = 0.24 and w e g — — 1. 

Given -ff(lna) and i?(lna), the metric ratio at super- 
horizon scales comes from solving Eqn. ([7]) 
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$ = 0, (fcjr-0). 
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We have used the f{R) relation [24 



S + * = _.b£^, 
H k H ' 



(24) 



which when combined with (f — and Eqn. (|15|) gives 



* = 



1 -B 



(k H -» 0) 



(25) 



The solution of Eqn. J24]) together with J25]) yields the 
metric ratio 

$ + $ 

g(lna, fc g = 0) = ffsH(lna) = _ . (26) 

The density evolution function can be adequately de- 
scribed by noting that (' k%BQ and that B also con- 
trols the behavior of g. We take 



ft = c <9 



(27) 



with cq k, —1/3. 

For the transition to the quasi-static regime we take 
the interpolating function 



g(lna, k) 



ffSH + ,9Qs(c g fcg)" 8 
1 + {c g k H )«« 



(28) 



where <7qs = —1/3. We find that the evolution is well 
described by c g = 0.71-B 1 / 2 and n g = 2. We show an 
example of this fit in Fig. [31 

Finally, the effective Newton constant is rescaled by 
fa and the quasi-static transition takes place near the 
horizon scale 



fa = fit, cr = 1 



(29) 



In Fig. [4l we show how well the PPF parameterization 
reproduces the full f(R) metric evolution for scales that 
span the Compton wavelength transition in a w e s = — 1 
and Bo = 0.4 model. We have checked that a wide range 
of f(R) models including those of [3(| produce compara- 
ble matches with these parameter choices. 



C. DGP Model 

In the DGP braneworld model, the transition between 
two different behaviors for g occurs at the horizon scale. 
Above the horizon, the propagation of perturbations into 
the bulk requires solving the full 5D perturbation equa- 
tions [3l|. Fortunately, above the horizon scale the evo- 
lution is scale free and can be solved using the iterative 
scaling method of [23|. Well below the horizon the evo- 
lution reaches the quasi-static limit where the equations 
can be effectively closed on the brane [(| [25[ . We there- 
fore take a similar approach to the f{R) case of interpo- 
lating between these two well-defined regimes. 
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FIG. 5: Evolution and scale dependence of the metric ratio g 
in the DGP model compared with a PPF fit. Here Q m = 0.24 
and the PPF parameter c 9 = 0.4. 



On super-horizon scales, the iterative scaling solution 
is well described by the fitting function 



gSH(lna) = 



9 



8Hr c 



1 



1 



0.51 



Hr c - 1.08 



(30) 



where recall that r c is the crossover scale. In the DGP 
model C' is again related to g and so we take to be de- 
fined by Eqn. (p?T)) with « 0.4. The expansion history 
is given by 



H 



where 
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(31) 



(32) 



For illustrative purposes we take O r 
In the quasi-static regime [25 1 



0.24. 



ffQs(hia) 



1 H' 

1 - 2Hr c ( 1 + 

1 3 H 



(33) 



We employ the interpolation function (|28[) to join the two 
regimes. In Fig. [51 we show a fit to the results of [23j with 
c g = 0.4 and n g = 3 for several values of ku that span 
the transition. The remaining parameters are 



/g = 0, 



Cr = l. 



(34) 



Around horizon crossing, the scaling assumption of [231 ] 
is briefly violated leading to possible numerical transients 
in the solution and an ambiguity in the exact value of c g . 
The results for the metric evolution in [23[ are best fit 
with c g = 0.14 as shown in Fig. [6) The transition pa- 
rameter c g should therefore be taken as a free parameter 
in the range c g ~ 0.1 — 1 until a more precise solution is 
obtained. 
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FIG. 6: Evolution and scale dependence of $_ in the DGP 
models compared a PPF fit. Here Q, m — 0.24 and the PPF 
parameter c g = 0.14. 



IV. NON-LINEAR PARAMETERIZATION 



FIG. 7: PPF non-linear power spectrum ansatz for an f(R) 
model. The non-linear power spectrum is constrained to lie 
between two extremes: defined by halo-model mass functions 
with the quasi-static growth rate [c n i = or Pa(k)] and the 
smooth dark energy growth rate with the same expansion 
history [c n i = oo or Poo{k)]. Here Bo = 0.001, w g = —1 and 
f2m = 0.24 with other parameters given in the text. 



As discussed in mi C[ we expect that a successful modi- 
fication of gravity will have a non-linear mechanism that 
suppresses modifications within dark matter halos. In 
this section, we construct a non-linear PPF framework 
based on the halo model of non-linear clustering. Al- 
though a complete parameterized description of modified 
gravity in the non-linear regime is beyond the scope of 
this work, the halo model framework allows us to incor- 
porate the main qualitative features expected in these 
models. Searching for these qualitative features can act 
as a first step for cosmological tests of gravity in the non- 
linear regime. 

Under the halo model, the non-linear matter power 
spectrum is composed of two pieces (see 32] for details 
and a review). One piece involves the correlations be- 
tween dark matter halos. As in general relativity, the 
interactions between halos should be well described by 
linear theory. The other piece involves the correlations 
within dark matter haloes. It is this term that we mainly 
seek to parameterize. 

Specifically given a linear power spectrum of density 
fluctuations Pl, the halo model defines the non-linear 
spectrum as the sum of the one and two halo pieces 



P(k) = Ji(fc) + ll(k)P L (k) , 



(35) 



with 
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where p — p m {\na — 0). Here the integrals are over the 
mass M of dark matter halos and dn/dhiM is the mass 
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FIG. 8: Fractional difference in P(k) of the PPF non-linear 
f(R) ansatz from the smooth dark energy prediction with 
the same expansion history. As c n i — » oo deviations become 
confined to the weakly non-linear to linear regime. The model 
parameters are the same as in Fig. [7] 



function which describes the comoving number density 
of haloes. y(M, k) is the Fourier transform of the halo 
density profile normalized to y(M, 0) = 1 and b(M) is 
the halo bias. Note that Ii(k = 0) = 1 so that the linear 
power spectrum is recovered on scales that are larger than 
the extent of the halos. 

A simple ansatz that restores general relativity in the 
non-linear regime is that the mass function and halo pro- 
files remain unchanged from general relativity. Specifi- 
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cally, whereas the mass function and halo profiles usually 
depends in a universal manner on a(M) the rms of the 
linear density field smoothed on a scale that encloses the 
mass M at the background density, we replace this with 
the rms of the linear density field of a smooth dark energy 
model with the same expansion history ctgr(M). For def- 
initeness, we adopt the Sheth-Torman mass function and 
bias [331 



dn 



dv 



= — f ( ) 



b(M) 



where v = 6 c /a(M) and 



1 



2p 



8 c [l + (av 2 )P] ' 



(37) 



vf(v) = AyJ^av 2 [l + {av 2 )- p ] exp[-ai/ 2 /2] ■ (38) 

We choose S c = 1.68, a = 0.75, p = 0.3, and A such that 
J dvf(v) = 1. For the halo profiles we take the Navarro, 
Frenk and White (NFW) profile 
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FIG. 9: Fractional difference in P(k) of the PPF non-linear 
DGP ansatz from the smooth dark energy prediction. The 
DGP model has less power than the equivalent dark energy 
model in the linear and weakly non-linear regime. Here Q m = 
0.24 and other parameters are given in the text. 



p oc 



cr/r vir (l + cr/r vir ) 2 



where r v ; r is the virial radius, I3E 



c{M v ) 



9 ( M 



l + z V M* 



-0.13 



(39) 



(40) 



and M* is defined as cr(M*) = S c . We call the result of 
taking a(M) = acji(M) in the halo model equations ([35]) 
Poo(k). 

At the opposite extreme, we can make the ansatz that 
the usual mapping of the linear to nonlinear power found 
under general relativity remains unchanged in modified 
gravity. In this case, changes in the linear growth rate 
determine the non-linear power spectrum. This type of 
prescription of adopting the linear to nonlinear scaling 
of general relativity has been tested against cosmological 
simulations with various modified Poisson prescriptions 
(36| . It represents the case where gravity is modified 
down to the smallest cosmological scales. Specifically, in 
our halo model we take a(M) = (Tppp(M) as calculated 
from the linear power spectrum of the modified gravity 
model and employ the same Sheth-Torman and NFW 
prescriptions as before. Let us call the power spectrum 
in this limit Po(k). 

We can parameterize an interpolation between these 
two extreme behaviors 



P(k) 



l + c nl E 2 (fc) 



(41) 



that is based on the degree of non-linearity defined by 
E 2 (fc)^^#. (42) 



2vr 2 



The analogous interpolation can also be used for the 
power spectrum of that enters into gravitational lens- 
ing observables [37| . 

We show an example of this non-linear ansatz for an 
f(R) model with fl m = 0.24, fl m h 2 = 0.128, VL b h 2 = 
0.0223, n s — 0.958 and initial curvature fluctuation of 
S c = A\, /2 = 4.52 x 10~ 5 at k = 0.05 Mpc" 1 . We fur- 
thermore fix the expansion history with w c ff=-i- 

Unfortunately no cosmological simulations exist for 
f(R) or DGP models against which to test the accuracy 
of this non-linear ansatz. Moreover even under general 
relativity, the halo model model of Eqn. (j35|) does not 
exactly reproduce the non-linear spectra of cosmological 
simulations. 

More robust in our parameterization is the relative 
change between a PPF power spectrum and the general 
relativistic prediction with smooth dark energy Pgr(&) 
and the same expansion history. This factor can then 
be applied to more exact results from cosmological sim- 
ulations to search for deviations from general relativity. 
The difference between Pgr(^) and P<x>(k) is that the for- 
mer uses the general relativistic linear power spectrum in 
Eqn. (|35p . This prescription can be further refined by cal- 
ibrating Po(k) directly from simulations of the modified 
Poisson equation [3(| . 

We show the fractional change between the PPF power 
spectra and the general relativistic power spectra for 
the f(R) model in Fig. [5] and a DGP model with the 
same parameters but with the DGP expansion history of 
Eqn. (|3"Tj) in Fig. [9l Note that as c n i — * oo, deviations 
appear mainly in the linear to weakly non-linear regime. 
For the f(R) model they appear as an enhancement of 
power and for the DGP model as a deficit of power re- 
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fleeting the opposite sign of g in the linear regime of the 
two models. 



V. DISCUSSION 

We have introduced a parameterized framework for 
considering scalar modifications to gravity that acceler- 
ate the expansion without dark energy. This framework 
features compatibility in the evolution of structure with a 
background expansion history on large scales, a modifica- 
tion of the Poisson equation on intermediate scales, and 
a return to general relativity within collapsed dark mat- 
ter halos. This return to general relativity is required of 
models to pass stringent local tests of gravity. We have 
also clarified the formal relationship between modified 
gravity and dark energy in the Appendix. A metric based 
modified gravity model can always be cast in terms of a 
dark energy component with a stress energy tensor de- 
fined to match its influence on the metric. However such 
a component would possess dynamics which are coupled 
to the matter. 

Our parameterized post-Friedmann framework fea- 
tures several free functions even in the linear regime. 
The most important function is the relationship g — 
($ + $)/($— "J*) between the time-time and space-space 
pieces of the metric in Newtonian gauge. Supplementing 
these are two functions that link the metric to matter 
density perturbations: one on super-horizon scales and 
one on intermediate scales. Finally there is a param- 
eter that controls the interpolation between these two 
regimes. 

We have shown that with an appropriate choice of pa- 
rameters this framework describes linear perturbations 
in the f(R) modified action and DGP braneworld grav- 
ity models. It may be used in place of the more com- 
plicated 4th order and higher dimensional dynamics ex- 
hibited in these models respectively when studying phe- 
nomena such as the integrated Sachs- Wolfe effect in the 
CMB, large-scale gravitational lensing and galaxy clus- 
tering. We intend to explore these applications in a fu- 
ture work. 

On non-linear scales our framework features an ansatz 
based on the requirement that scalar modifications 
should be suppressed locally in order to pass the stringent 
tests of general relativity in the solar system. Indeed the 
scalar degrees of freedom in both the f(R) and the DGP 
models possess non-linearities that drive the dynamics 
back to general relativity in high curvature or high den- 
sity regimes. Our ansatz is based on the halo model of 
non-linear clustering. It allows for a density dependent 
interpolation for the abundance and structure of dark 
matter halos between the expectations of general relativ- 
ity and the modified Poisson equation on intermediate 
scales. 

Due to the current lack of cosmological simulations in 
these modified gravity models, the accuracy of our simple 
ansatz remains untested. With cosmological simulations, 



our framework can be extended and refined by introduc- 
ing more parameters that describe the potentially mass- 
dependent modification of dark matter haloes. In fact, 
our simple halo model parameterization is not even suffi- 
cient to accurately model non-linear effects in general rel- 
ativity. Nonetheless phrased as a simple template form 
for relative deviations in the power spectrum between 
modified gravity and general relativity with smooth dark 
energy, our current ansatz can be used in conjunction 
with more accurate results from dark energy cosmologi- 
cal simulations. For example, it can be used to search for 
possible deviations of this type as a consistency check on 
dark energy inferences from expansion history tests with 
upcoming cosmic shear surveys. 

While many such consistency tests have been proposed 
in the literature, it is important to incorporate a density 
dependence to the modifications as we have done here. 
The principle that non-linear scales should exhibit a re- 
turn to general relativity itself suggests that mildly non- 
linear scales provide the most fruitful window for cos- 
mological tests of gravity. Furthermore uncertainties in 
the baryonic influence on the internal structure of dark 
matter halos in the de eply non-linear regime even under 
general relativity (e.g. [38|,[39j]) make consistency tests in 
this regime potentially ambiguous. Our parameterized 
framework should enable studies of such issues in the fu- 
ture. 



APPENDIX A: DARK ENERGY 
CORRESPONDENCE 

Suppose we view the modifications to gravity in terms 
of an additional "dark energy" stress tensor. We are free 
to define the dark energy stress tensor to be 



Given this association, all of the familiar structure of 
cosmological perturbation theory in general relativity ap- 
plies. In particular, covariant conservation of the matter 
stress energy tensor Tf£ and the Bianchi identities imply 
conservation of the effective dark energy 

V M TT = . (A2) 

The remaining degrees of freedom in the effective dark 
energy stress tensor can then be parameterized in the 
same manner as a general dark energy component [H, [4l| . 
Two models that imply the same stress energy tensor of 
the combined dark matter and dark energy at all points in 
spacetime are formally indistinguishable gravitationally 

Note however that this equivalence is only formal and 
two physically distinct models, e.g. f(R) modified grav- 
ity and scalar field dark energy, will not in general imply 
the same effective stress energy tensor [43[ . The Einstein 
and conservation equations do not form a closed system 
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and the distinction between modified gravity and dark 
energy lies in the closure relation. For dark energy that 
is not coupled to matter, the closure relationship takes 
the form of equations of state that define its internal dy- 
namics. These micro-physical relations do not depend 
explicitly on the matter. For example for scalar field 
dark energy, the sound speed or the relationship between 
the pressure and energy density fluctuations is defined in 
the constant field gauge without reference to the matter 
[H , and is associated with the form of the kinetic term in 
the Lagrangian [44| . 

For modified gravity of the type described in this pa- 
per, we shall see that the closure relations must depend 
explicitly on the matter (see also [HJ). The effective dark 
energy of a modified gravity model must be coupled to 
the matter. In other words, while the modification of 
gravity can be modeled as fifth forces mediated by the 
effective dark energy, it cannot be viewed as a missing 
energy component that obeys separate equations of mo- 
tion. 

It is nonetheless useful to phrase the PPF parameter- 
ization in terms of an effective dark energy component. 
It enables the use of the extensive tools developed for 
cosmological perturbation theory and facilitates the de- 
velopment of PPF formalisms in different gauges. 



1. Covariant Field and Conservation Equations 

Following (2lL [45I ] , we parameterize linear scalar metric 
fluctuations of a comoving wavenumber k as 



n a J ,l 2.1V 1. 

l~ 

a- 2 ( 7 y - 2H L Yj ij - 2H T Y ij ) , 



g Qi = -a- 2 BY 1 , 



.9 



(A3) 



where the "0" component denotes conformal time 77 = 
J dt/a and 7^ is the background spatial metric which we 
assume to be flat across scales comparable to the wave- 
length. Under this assumption, the spatial harmonics are 
simply plane waves 



Y = e ik x , 

Yi = (-fc)v l r, 

-2 



Yn = (AT 



ViV, +~/ ij /S)Y. 



(A4) 



matter plus effective dark energy stress tensor. For sim- 
plicity we assume that the radiation is negligible during 
the epochs of interest. 

By definition, Eqn. (|A1|1 enforces the usual 4 Einstein 
field equations [22|] 



1 



B 

3"' ' k H 
AttG 



H f rp 



H 2 k 2 H 



5p + 3(p + P )- 



B 



A + Hl H — -— 



3~ 

87TG 

lpk% 



B' + 2B 


HtL 




H' T 


k H 


. K H 




k 2 



pU, 



. TTl EL 4ttG. .v-B 



A' 



' H 



HI 



k 2 
3 

W 
V H 



A-^L(B' + B) 



H' L = 



AttG 



1 



(Sp + fSp), (A6) 



where recall ' = d/dlna and kn = (k/aH). The con- 
servation laws for the matter and effective dark energy 
become 



8p' + 3(5p + 5p) = -(p + P )(k H v + 3H' L ) 



(A7) 



{aHp + p)(v-B)}' 
a 4 k H 



= Sp- -pU + (p 



p)A. 



There are 4 metric variables and 4 matter variables per 
component that obey 4 Einstein equations and 2 conser- 
vation equations per component. However 2 out of 4 of 
the Einstein equations are redundant since the Bianchi 
identities are automatically satisfied given a metric. Fur- 
thermore, 2 degrees of freedom simply represent gauge or 
coordinate freedom. This leaves 2 degrees of freedom per 
component to be specified. Usually, this involves defin- 
ing equations of state that specify the spatial stresses in 
terms of the energy density and velocities. As we shall 
see, it is this prescription that must be altered to describe 
modified gravity. 



Likewise the components of the stress tensors can be pa- 
rameterized as 



2. Gauge 



T\ = -p-5p, 



T l = -(p + p)vY\ 

T) = (p + 5pY)S l J + P UY) 



(A5) 



The scalar gauge degrees of freedom are fixed by gauge 
conditions. Under a gauge transformation defined by the 
change in conformal time slicing T and spatial coordi- 
nates L 



where we will use the subscripts m to denote the matter 
and e to denote the effective dark energy. When no sub- 
script is specified we mean the components of the total or 



V = V + T , 
x* = x* + LY l . 



(A8) 
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the metric variables transform as 

A = A - aH(T' + T) , 
B = B + aH(L' + k H T) , 

H L = H L -aH(T+^k H ), 
H T = Ht + aHk H L , (A9) 
and the matter variables transform as 

Sp = Sp — p'aHT , 
Sp = Sp — p'aHT , 



v = v + aHL' , 

n = n. 



(A10) 



A gauge is fully specified if the functions T and L are 
uniquely defined. 

In this paper we work in the matter comoving and 
Newtonian gauges. The matter comoving gauge is spec- 
ified by the conditions 



B = v m , 
H T = 0. 



(All) 



They fully specify the gauge transformation from an al- 
ternate gauge choice 



T = (v m - B)/k , 
L = -H T /k . 



(A12) 



To avoid confusion between fluctuations defined in differ- 
ent gauges, we will define 

e = a, 

P A = Sp, 
Ap = Sp , 
V = v. (A13) 

Ap should not be confused with pA = p(Sp/ p). 

The appropriate Einstein and conservation equations 
for this gauge can be obtained by utilizing these defi- 
nitions in Eqn. (|A6[) and (|A8[) . For example, the third 
Einstein equation reads 



C = — —\Pz+Pz)- 



H 2 



(A14) 



and the energy-momentum conservation equations for 
the matter become 

A' m = -k H V m - 3C , 

£ = 0. (A15) 

The dark energy momentum conservation equation 

[a*(pe+Pe)(Ve-V m )]' * 2 

-5^ = A Pe - - Pe U e , (A16) 



in conjunction with Eqn. (|A14[) and the first Einstein 
equation implies that unless Ap e or p e Tl e > 0(Ap/kf 1 ), 
C'/C -»■ as k H -»• 0. 

Similarly, the Newtonian gauge is defined by the con- 
dition B = Ht = and the transformation 



kk 



H 



L = — 



To avoid confusion we define 



(A17) 



A. 



(A18) 



We refrain from utilizing matter variables in Newtonian 
gauge but note that velocities in the two gauges are the 
same. The relationship between the two metric fluctua- 
tions are 



k H 
V 



(A19) 



The matter momentum conservation law in Newtonian 
gauge becomes 



vL + v m 



k H m . 



(A20) 



This equation can alternately be derived from the gauge 
transformation equation (1A19|) given that £ = 0. 

Finally the Einstein equations (|A6|) in Newtonian 
gauge imply 



$ + * 
2 

$ - * 



WkT pn 



H 2 B Pe e ' 



rr [ pA + Z{p + p) 



H 2 k 



H 



AttG 



p m A m + p e A e + 3(p e + p e ) 



pli 

V, 



(A21) 



Vm 



(A22) 



where we have assumed that the anisotropic stress of the 
matter is negligible. A finite metric ratio parameter g = 
+ — vf) is thus associated with a non- vanishing 

effective anisotropic stress. 



3. PPF correspondence 

The system of equations defined by the field equations 
and the conservation equations are incomplete. To close 
the system of equations two more conditions must be re- 
quired of the effective dark energy. It is this closure con- 
dition that the PPF parameterization must determine. 
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Given that the matter has no anisotropic stress, 
Eqn. (|A21[) defines the anisotropic stress of the effective 
dark energy in terms of the metric 



AttG 



(A23) 



where recall $_ = (<!> — ^)/2. This is the first of two 
closure relations. 

The second closure relation comes from equating 
Eqn. (|A22[) and the modified Poisson equation fTf 



p e A e + 3(p e +p e ) Ve — +p e II e 



Cff 



47rGa 2 



r. (A24) 



The PPF equation of motion (|2lj) for T is therefore the 
"equation of state" for the effective dark energy. 

The conservation laws for the effective dark energy 
and/or remaining Einstein equations then define the 
other two components V e and Ap e . For example, 



V e = Vrr, 



H 2 



4i:Ga 2 (p e +p e ) 



(A25) 



1 



Ap e = Pe A e - - Pe A' e - (p e + Pe )(k H V e /3 + C') ■ 



The modification represented by this prescription obeys 
all 4 Einstein equations and both sets of conservation 
laws. 

Unlike the case of a micro-physical candidate for dark 
energy such as a scalar field, the closure relations not 
only cannot be defined as direct relationships between 
the spatial stresses and the energy density and velocity, 
they here involve the matter and the metric fluctuations 
directly. Hence the effective dark energy is implicitly 
coupled to the matter and cannot be described as an 
independent entity. 

On the other hand, the virtue of making this cor- 
respondence explicit is that with these relations all of 
the usual representations of perturbation theory can be 
reached by standard gauge transformations from our 
matter comoving and Newtonian representations. 
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